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SUMMARY 



Three general theories treating the potential flow 
about an arbitrary wing section are discussed in this re- 
port. The first theory treats the method of conformal 
transformation as laid down by Theodorsen and Garrick; the 
second is a generalization of Birnbaum's theory for mod- 
erately thick airfoils; the third is a general investiga- 
tion of the complex velocity function with particular ref- 
erence to the relations first discussed by F. Teinig. 

The relative merits of the different methods in ques- 
tion are illustrated on a worked-out example and will be 
published in a subsequent issue of this periodical. 



INTRODUCTION 



The present investigations relate to the two-dimen- 
sional potential flow of a f rictionless , incompressible 
fluid around any simply connected region, particularly 
around airfoils. Such a flow, as is known, is completely 
de scribed by a regular, analytical function - the complex 
velocity function. Since an analytical function, in turn, 
is completely described by its values along any closed 
curve I remaining wholly within the region of regularity, 
it is sufficient to know this complex velocity function 
along a curve enveloping, once, the zone washed by the 
stream - for which purpose the limiting curve or profile 
contour itself may be chosen, with observation, of course, 
of any potential singularities on the limiting curve. In 



"Theorie der ebenen Potential stromung urn beliebige Trag 
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any case the pro Diem may "be considered solved, once the 
velocity along the profile contour is known in magnitude 
and direction. 

The very next problem is to find the flow - that is, 
a"bove all else, the velocity along the contour for any- 
predetermined region. In this case the contour with its 
direction affords the direction of the velocity as well, 
leaving simply the quest of the function for the absolute 
value of the velocity along the contour. This, as is known 
(reference l), requires the supplementary assumption of: 

1) the velocity at infinity in magnitude and direction 4.1 otu\ 

2) either the circulation, i.e., the line integral of the 
velocity along a curve encircling the obstacle once, or 
else the point on the contour of one of the stagnation 
point $3 of the flow. 

This general flow problem can be reduced in known man- 
ner to the mathematical problem of conformal transforma- 
tion of the contour of the washed region onto the contour 
of a circle - a problem whose solution is afforded by 
Riemann's law of transformation. It can be considered 
solved if it is possible to construct this conformal func- 
tion not merely by visualization but by actual plotting, 
according to a method which must be rapidly converging, if 
infinitely many steps are necessary. Such a method was 
advanced and proved by T. Theodorsen and I. E. Garrick 
(reference 2) and in the treated cases, yields very quick 
results. It can also be proved that this method converges 
under certain sufficing assumptions and solves the trans- 
formation problem rigorously. The results, built up on 
it, are the mathematically exact solutions of the general 
flow problem. 

But in aerodynamics, the opposite of the above problem 
is also of significance - i.e., to find a profile along 
which the flow .about the contour is accompanied by a de- 
sired velocity distribution. ■ F . Weinig has established 
the surprising fact that this- problem is mathematically 
simpler* and that, if the desired, velocity distribution in 
the potential flow about a simply connected region can 
take place at all, the contour of this region can be ex- ■ 
actly determined without an infinite method. -Admittedly, 
the assumptions required for the velocity distribution 
along the contour itself cannot be summarily complied with, 
hence it is impossible to' take such distributions directly 
as a starting point for the calculation of profiles. A 
suitable assumption, instead of this velocity distribution, 
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is a freely available real function following from it by 
'distortion in abscissa ^direction. The method is particu- 
larly suitable for- .developing profiled with prescribed 
characteristics. 

In. earlier airfoil theories', special conformal funo- 
tions with a .varying number of parameters which give the 
conformal transformations "between circles and airfoil- 
shaped regions, formed the center of investigation, and 
logically the results were predictions about special, multi- 
parameter airfoil families. General relations "between the 
functions for the profile contour and for the velocity dis- 
tribution could not he obtained in this manner, with the 
notable exception of infinitely thin airfoils by the so- 
called "Birnbaum the o'ry " ' whi ch , proceeding from the image 
of a nonuniformly covered flat vortex layer, approximately 
established a relation between the difference in velocity 
on the upper and lower profile surface and the slope of 
the mean line of the profile. Notwithstanding its re- 
stricted range of validity, this theory presented many ad- 
vantages and made its generalization to include thick air- 
foils very desirable which, however, did not succeed sat- 
isfactorily with the vortex layer concept. 

In the present report it is attempted to combine the 
past information on potential flow around airfoils, so 
far as they are of this functional type, with a view to 
obtaining data which tie the profile contour and the com- 
plex velocity, or velocity distribution along the contour 
in plainest and most amenable form for calculation. 

In this connection, the findings by Theodorsen and 
Weinig are significant. It might "be surmised that by van- 
ishing profile thickness these theories lead back to the 
old theory for thin airfoils. *But such is not the case. 
It rather affords a new approximate theory - probably 'suf- 
ficient for all cases encountered in practice which, with- 
out restricting itself to infinitely thin profiles, con- 
tains the equation part of the theory of lifting vortex as 
limiting case and to that extent represents its general- 
ization. It ultimately affords a survey over a system 
of integral equations, every one of which is involved as 
mathematical starting point for a. general airfoil theory. 
It is shown how the theory pf Ixfting vortex surface and 
Weinig's results aline- themselves in this general arrange- 
ment'. 
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1. THE POTENT IAL ■ -MOW: AROUND A C IRCULAR CYL INDER AND 
ITS RELATION TO- FLOW AROUND GENERAL PROFILES 



The physical conditions of two-'dlmeh s ional vortex- 
free flow in frictionless, incompressihlfe. fluid are con- 
tained in the two equations ' 

^ = |u + |2C ■= 0 and |3 - |£ = 0 . (l ) 

^ dx 3 y dy dx . . • 

where u and v indicate the velocity" in ' x and y 
direction. They represent the Cauchy-Riemann differential 
equations for the analytical function' \ 

jU, . w*(z) = u(x,y) - 1 v(x,y) = (2) 

7% " . . a * ' ^ 

the so-called "complex ye lo.city 'function" which completely 
describes a two-dimensional potential flow. 

If it relates to the flow around a finite region, as, 
for instance, a circle or airfoil in unlimited fluid, sev- 
eral general predictions as to function w (z) can "be 
made. Since infinitely great velocities can ocgur only on 
the "border of the fluid region, the function w (z) is, 
ahove all, regular everywhere outside the region. Hence 
for grljft z, w*(z) can he represented hy the series 



w*(z) = O 0 + -i + % + ... (2a) 
z z 



The constant C Q indicates the complex velocity at infin- 
ity.. The 'constant Ci is purely imaginary, if the con- 
tour around the ohstacle is a streamline, as presumed 
here, in ,accord with classical theory. 6 1 is associated 
.with the circulation V of the "body through the relation 
r 

C x = i --• While ' C-, is decisive for the force exerted 
2tt 

hy the flow on the "body, C 2 is decisive for the moment 
on the "body." If C n =' A^ + i B n and p is the density 

of the flowing medium, .fframmel's well-known equations (ref- 
erence 3) give for the two components P^. and Py of the 
force and for the moment M turning clockwise about z=0: 
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P x .= 2iT.p B 0 B lt P y =' 2tt p A 0 B 1) - M = - 2tt p(A 0 B 3 + B 0 A a ) 

' : ' • ' • (3) 

For the most general flow about the unit circle 

with the two stagnation points z x = e *(a-0) and z ^ _ 

giCTr+a+P)^ tlxe com pi e:!C velocity- function is: 

w (z) = ¥ e ~ la 1.1 - S — . 1 + ---- — ~ (4) 

The flow has the ah solute velocity W at infinity where 
it forms the angle a with the x axis. Transformed, 
equation (4) gives: 

w*(«> = * (e-ia + 2i_3in_3 _ e^ 

The series (2a) stops in this case with the term z~ 2 . 

The coefficient of z" 1 is, as should "be , purely imaginary. 

On the circle circumference z = e iep it is 
w*(cp) = 2i W e - icp (sin (cp-a) + sin 3) 
hence the absolute velocity on the circle periphery is: 

|w| = |w*(cp)| = 2W | sin (cp-a) + sin B| (5) 

The velocity is zero at the two stagnation points for 

cp = a - 3 and cp ' = tt + a + 3. For the direction of flow 

at point z = e icp , we find u/v = - cot cp, which confirms 
the fact that the function (4) actually represents a poten- 
tial flow about the unit circle. 

The potential flow about the circular cylinder is, as 
known, of fundamental importance for the general theory of 
potential flow about any simply connected region B, be- 
cause the conformal function "z(z) which transfers the 
contour of the unit circle on the contour of B, makes it 
possible to deduce the velocity function w .("z^ for a flow 
around region B from the complex velocity w (z), accord- 
ing to equation (4). The relation is: . 

w*(z) = w*(z) A* ■ (6) 
dz 
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The- derivation of equation (.6'.) usually proceeds from the 
fact that in the conformal transformation by "z(z) the 
orthogonal fields of the streamlines and the potential 
lines of "both flow.s merge. ffinJliv*' e»-lr«^v c. -* iMtaaX^i 

JiwUk, .. . t -> . ' " i. 

But this formula can also he understood without re- 
sorting to the flow potential". For w ("z ) is an analyti- 
cal function, for whose functional values along the "bor- 
der of B the negative argument agrees with the direction 
of the contour. Function w'(z) meets this requirement 
for the unit circle. In the conformal transformation of 
"z(z), point z moves toward "z and the vicinity of z is 
turned through an angle ^arc d"z/dz. Hence, to insure the 
required direction. of w (z) at the horder of B, it is 
nece s sary that : 

- arc w*(z) = - arc w (z) + arc -f- 

dz 

But the analytical funcfirn which complies with this re- 
quirement is w*(z) 

dz 

This pro.cess affords -the potential flow- about an infi- 
nitely thin, flat plate. The analytical function which 
transfers the outside .of the unit circle in' the z plane 
into the ~z plane rectilinear from -1 to +1 reads: 

+ or z = ¥ (t)/**- 1 (?) 

which, written for z in equation (4). and multiplied "by 
d.z/d.7,, according to equation (6), gives for the complex 
velocity of the most elementary flow a"bout the flat plate 
extending from' -1 to +1' the following (with abbrevia- 
tion f = 217) : 

w*(z) =1. (cos a. i sin a - --1 . + ■ l ^LS-iJ ) ( 8 ) 

^ SI* - 1 - si^n J 

■ This equation describes a three-parame ter_sy stem of 
potential . flows . ■ Two v that is, the velocity W at infin- 
ity and angle of flow a can be- regulated at will.-. The 
value of the third -parameter ■ - 3 is physically conditioned - 
i.e., the Joukow ski ; condit ion of finite velocity at the 
trailing edge z" = 1 . In the present case w (l) is finite 

if lim (z sin cc - sin. 3).= 0, that is, if 3 = a. Thus, 
z >1 
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the velocity function for the flow around the flat plate 
f inally reads: - . 



w*(z) = . W (cos a - i sin a J =~~j) 



(9) 



In order to describe the velocity distribution on the 

surface of the, plate, the variable cp is employed. Ac- 
cording to equation (-7) "2 = cos cp = £" for * = e it; P. Then 
it is for the surface of the plate: 

w (cos cp)=W (cos a+sin a J ) = ^ (cos a+sin a tan |Q - 

after which the absolute velocity at point x = cos cp "be- 
comes 

|w(cp)| = 1 | cos a + sin a tan |*j (10) 

With this method any number of other potential flows 
can be mathematically described by different, choice of 
confo.rmal function. As it is easy to give conformal trans- 
formations which transform the contour of the unit circle 
into that of a simply connected airfoil-like region, it is 
equally possible to give profiles which may be controlled 
by mathematical theory. 

The theories of general potential flow, about airfoils 
discussed hereinafter, give the relation between the func- 
tion for the profile contour and the complex velocity along 
the contour, and tie these functions in a fashion amenable 
to calculation. The start is made with the strict formu- 
las for the flow about any predetermined airfoil to which 
the theory of Theodorsen and. Garrick leads.- 



2. RIGOROUS THEORY OF POTENTIAL FLOW ABOUT ANY PROFILE 



The previously described classical method of obtain- 
ing the potential flow about any simply connected region 
yields - in conjunction with the method of Theodorsen and 
Garrick f«r obtaining the conformal function for any ini- 
tial region - a rigorous solution of the basic problem of 
airfoil theory; that is, to find the potential flow about 
any given airfoil (references 2 and 4). 
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If the values of the complex velocity along the pro- 
file contour are known the pro"blem is solved-, "because 
then Cauchy's integral formula gives the velocity func- 
tion for any point on the outside region replete with the 
flowing medium. To 'find' the velocity distribution on the 
profile contour, the procedure is as follows: The con- 
formal transformation (7) which transfers the slotted 
plane to the outside zone of the unit circle, transforms 
the profile into an almost • circular region. With the aid 
of the function which conformally transforms the contour 
of thi.s region onto the contour of the unit circle , the 
velocity distribution on the contour is then computed for 
the flow around this region. Prom this distribution -the 
velocity distribution along the profile contour is deduced 
with the aid of the first conformal transformation. 

Assume the arbitrarily given profile in the ~i plane 
to be so plotted that its trailing edge coincides with 
point z = 1 and point z = - 1 lies in the profile near 
the nose of the profile (fig, l). If the " z plane is con- 
formally transformed by function ' £ = "z .+ J z a - 1 on the 
t plane, then for the given position of the profile with 
respect to the fixed points "z = ± 1 , the image curve • of 
the profile contour in the I' plane is the contour of a 
simply connected region which; for the common profile 
forms, is in more or less satisfactory agreement with the 
unit circle. 



periodic with the period 2tt;' the boundary curve, viewed 
from the origin £ = 0, is star-shaped. 

The problem is to'' transform the contour of the unit 
circle z = e ic P on the contour of the region bounded by 
i(Q). This conformal transformation is standardized by 
the condition that the infinitely remote point of the z 
plane and the direction of the positive real axis in it, 
are coordinated to infinitely remote point of the £ plane 
and the direction of the positive real axis. According to 
Riemann, there is exactly one analytical function £(z) 
which meets this requirement and conformably transforms the 
two outside regions. But in view' of this fact, the con- 
formal function £ ( z ) has the form. 



The boundary curve 
(fig.. 2). The function 




I (z ) = z e 



F(z) 



N.A.C.A. Technical Memorandum No. 886 



9 



z 



where F(z) for |z| > 1 inclusive of z = <», is a reg- 
ular analytical function which assumes a real value for 

Z *. CO. 

Owing to z = e ic P and £ = e^ +i6 , it is 

on the unit circle. Thus the function F(z) assumes on 

the circle z = e iC(:> the boundary values \\f + i(8 - cp ) . 
Given either 8 (cp ) or ^(cp), the function F(z), and 
hence the conformal function £(z) can "be directly writ- 
ten. 1 Since F(z) takes a real value at infinity, we have 
here h(csa) = 0; hence "by steady ^ (cp ) for any z with 
|z| > 1. the equation 



If F(z) ss g + i h is an analytical functign that is reg- 
ular for |z| > R including z = m and -frCe-affy on the cir- 
cle |z| = E, all z = i? e icp with u > E follow the Pois- 
son integral formulas: 
2tt 

g<v,cp) = f g(a,cp«) Jff-r-Sf — _ _ aqp', 



2tt 



v - 2u E cos(cp'-cp) + E c 

2TT 



g (oo) = J- y° g(E.cp') dcp' 



and 2tt 

g(u,cp) = g(oo)- JL / h(R,cp«) SV-I-Sinlcej^cpi dcpI 

,£Tr '-' u 3 - Zx> R cos(cp'-cp) + E 

o 

which, combined, read: „ 

n icp» 



P(z) = i h(.) + i / g(E.cp') dcp t 



The second equation yields for E = 1 , v — >E, the impor- 
tant formula 2tt 

g(l,cp) = g(oo) - JL h(l,cp' ) cot ^-f-^ dcp' 

The £-irg%T*%4Jste integral here stands for Cauchy's principal 
value (cf. Harry Schmidt, Aerrdynamik de s Tluges, p. 87). 
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2tt 

T(z)=± /\<q>> ii4j d<P (ID 
./ z — e ^ 

o 

is applicable, while on the unit circle itself, we find: 

2tt 

- cp + G.(cp) = ~. ^ t(tp') cot - — |— dcp' 

" ' '" '"' 0 

Introducing the phase difference cp - 8(cp) as new 
function e(cp) 

. ' €(cp) = cp - 6(cp) • • .... (12) 

gives 'for the present 

2tt 

e(cp.) .= =1 I \|/(cp') cot dcp' 
.... 2tt . / -2 



\ 



If ^(cp) has a continuous first derivative, partial inte- 
gration leads to 

2tt 

. . e.'(cp) ■= - 1 - / in sin 2 5lL=-SP ' d cp ' 

2tt .y dcp 2 

o : -. ■■• •• 

Actually this is not the function $(q>) "but rather 
\i/(6). and the prohlem is to find the function e(0) for 
the given \|/(9). Then the function cp(9) follows from 
e(9) .according to equation qp(0 ) = G + c(6), and cp(6) 
together with ^(6) gives the function \|/(cp). With \j/(cp) 
the conformal function can "be explicitly written, accord- 
ing to equation (11 ). 

There remains the relation "between ^(6) and e(G) 
which, howeyj^r^J^s easily established when assuming that 
dcp/ d6 is STt-e-a-dy and other than zero, "because then 



V 
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2TT. 

J <lcp 2. ■• ^ 



2lT 



In sin 3 - 

d.0 



and consequently, 



c(e«) = JL 



2TT 



[e'-e e(G')-e(e)' 



d6 




sin' 



2 2 



de 



(13) 



This is the equation upon which Theodorsen's meth- 
od is based. It is a nonlinear and singular integral equa 
tion and therefore does not aline itself in the known the- 
ory. But, according to Theodorsen, it can he successfully 
applied by iteration process. This method is similar to 
the Picard-L'indelof method of solving differential equa- 
tions. It consists in temporarily substituting the unknown 
e(8) on the right-hand side for any suitable function 



(for instance, e 0 (6) = 0). Then equation 



2TT 



€ X (9) = 



JL 

2tt- 



/ M 

/ d8 



In sin e 



6J_- 
2 " 



de 



(13a) 



defines £i-(6) and generally gives e k+].(8) 
of e k (8 ) through the recursion formula: 



with the aid 




sin c 



[*LL=_£ + £ki6_0_r_£k(e_) 



de 



P-o-r— Jb-r-o-o-f- 



rrt— 4h-e 



convergence it would have to be 
shown that the functions e k (9 ) tend toward a' limiting 

function which satisfies the integral equation (13). This 
proof is not adduced in the original work of Theodorsen 
and Garrick, but rather in the practicability of the de- 
scribed method instead, by showing on several solved exam- 
ples that in them the functions € v (6) are practically no 
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longer di st ingui sha Die after a few iterations. The solu- 
tion of the integral is effected with a quadratic formula 
which considers 20 points in the interval (0^2 tt ) . 

Incidental to the present work, some searching inves- 
tigations concerning the - convergence of this method were 
carried out. It could "be proved that the functions e, (8) 

converge uniformly toward a &t-e«r4y limiting function e(0), 
which complies with equation (13), if the initial function 
\j/(6) can "be represented by a Fourier series whose coeffi- 
cients a n and b n are of the order of magnitude of 

A M- ■ 
, where p > 2 and A a^e smaller than an individ- 

(n+l) P 

ually specifiable, -positive figure - A Q different from 0. 
But the calculations of this proof are too extensive to-be 
reproduced here. ■ 

It might be noted that the very first step of the 
method affords, in general, such a satisfactory approxima- 
tion that it suffices for practical purposes. The numeri- 
cal evaluation of the integral of equation (13a) may be ef- 
fected with Theodorsen's calculation scheme. 

Suppose that the function e(8) is known: then, be- 
cause cp(0) = G + 6(6) and = 1 + |~g , according to 
equation (ll), the conformal function £(z) is: 



2u 1 W) ;:;TTe-+e7 k 1 + di; 



d6 



Uz) = z e 0 . * c (14) 

from which the enlargement ratio of the conformal trans- 
formation at infinity follows immediately at 

2tt 



a = lim = e 2Tr 9 



But on the boundary of the unit circle we have, because of 
t = e^ +ie at the point' of z = e 1 ®* 6 * = e ic P 



N.A.C.A. Technical Memorandum No. 886 13 

(§£\ - *e _ We * V e 1 -de *_ ie 
=W ; z=ei cp - e iXe + Tj - a -( x + H) e iTe+Tr i + |f 

(15) 

By reason of this the' Velocity distribution in the 
t plane is possible. The values of the complex velocity 
function w**(£) along the contour follow from w*(z), 
according to equation (4), and d£/dz according to 

w**(8) = «"(t) tw ^ l9 = (»•(.) ff) z=eW> 

But on the circle periphery it is, according to equation 
(4) 

w *(z) _ Q icp = 2i W e- ic P [sin(cp-a) + sin 6] = 

2i ¥ e - i ( e+£ ) [sin(£+e~cc) + sin B] 

Hence 

1 + ^1 

w**(6) =2il e-^-iQ |f [sin(6+c-a) + sin B] (16) 

1-i 

de 

For the determination of the complex velocity w ("z) 
in the "z plane along the profile contour, the transfer from 
plane £ to plane "z must "be effected. Both planes are mu- 
tually related through the conformal function" £ = ~z + 

<J z 3 - 1 or "z = ^ (^t + jr^. Using the ellip'tic coordi- 
nates \|/ , 0 the equation of the profile contour in plane 
z reads : 

z(0)=§ 0 + ^_^ +ie =|[ei(e-^) + e-i(e-ivl/).] „ CO s(0-it) 
Along the contour, it is: 

(&) = 1 (t 1\ = e:^r^r e i(0-ixl/L e -i(0-i\lr)n , 

«= i e-^- 3 - 6 sin(0-i\|/) 
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Besides the enlargement ratio of this transformation at 
infinity is ^ lim^ (jjj^) ~ \* hence, the velocity' distribu- 
tion along the' profile contour: 

-, . d e 

w*(e) =2W « ■ *5l£ +€-») + a in. P (17) 

- . d\|/ sxn(6-i\k) 

Now the parameter 3 must "be so defined that at the 
trailing edge, that is, for 6=0, the velocity remains 
finite. The condition for this is: 

3 = a - e 0 with e Q = e(Q) Q = a 

To make the final result agree with the result for the 
flat plate, equations (9) and (10), the insertion of this 
value for P is followed "by the following changes: 

sin(G + e -a) + sin(a~e 0 ) = sin(a-e 0 ) + sin [(6 +e -e 0 )- (a-e 0 )] = 
= sin (a- % ) + sin ( G+e~ e 0 ) co s (a- e 0 ) - co s (6 + e- e 0 ) sin(oc-e 0 ) = 

= sin(6+e~e n ) 



, 1-cos (6 + e-e 0 ) 
^cos(a-e n ) + fl in(o-e 0 ) ~J-J^Z1~Y J 

= sin(G + e-e 0 ) j cos(a-c Q ) + sin(a-£ 0 ) tan---™ 1 

In addition, we substitute the absolute flow velocity 
W in the plane i for W. It is with the constant a 
according to equation (14a) 



z 



lim 

giving as final profile contour 

i + 

w *(6) = a W — ------ L 

d6' 



w*(z) -4-11= W - and hence 2W = a¥ 
d£ dz | a 



c os ( a- e c ) + 



+ sin(a-e 0 ) tan ^IIzLsl | ( 18 ) 
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But, on account of 



|sin(e-I*)| = J sin 2 € + sinli \ and | 1-i || | = "j~ 



1 + 



the absolute velocity, takes the form 



|w(e)| = a I 



sin(6+g~e 0 ) 



+ d 6 



V 



1 + 



[ 



cos(a-« 0 ) + sin(a-€ 0 )tan ^ — °- j 



(19) 



Equations (18) and (19) are the strict equations for 
the velocity of the potential flow along the profile con- 
tour. To obtain an explicit expression for the complex 
velocity at any other point z of the outside region, the 
Gauchy integral formula is applied. Since the complex ve- 
locity function w*("z) outside of the profile inclusive 
of infinity is regular throughout, if "z denotes any point 
of the outside zone and "z f a point of the profile con- 



tour, the equation 2 reads 
W* ("i") = w(eo) + 



2-rri y C i - i' 



= W eT ia at infinity, and z » = cos (9-1^), 
hence d"z 1 = fl - i ^ sin(e-i x i / ) d6 on the profile con- 

\ d e s 

tour. Thus equation (17) "gives the complex velocity func- 



tion in the form 



2tt 



_ ._ iaW P sin(6+e-a) + sin(a- ^ ) 



If | z| > |cos(6-i\l0| , the series development for 
("z) is " . 



(20) 



If F(z) is a regular function of a simple region G 
"bounded at infinity and everywhere else outside of a part- 
ly smooth curve C, and T?(z) is' steady on' 0, every z 
of the outside zone of G follows Gauchy ' s integral for-. 



mula 



F(z ) = F(oo) + 



2tt i 



P( i) 



See also reference 1, p. 68. 
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C 0 

• w*(z). = W e~ ia + =- + =1 + • . • 1 '• 

with 2tt 

c n = |^ [sin(G+e-a> + sin (a- e 0 ) G] co s n " 1 (G -i v|/) 

1 ( 1 + H)« ' 



► (20a) 



This integral can "be computed for the coefficient G 1 and 
leads to a known formula. It is with Q + € = cp 

2TT 

Ci'= irf ^T[sin(G + e-a) + sin(a-e 0 )] (l + £§) dG = 

0 2TT 

= W / [sin(cp-a) + sin(a- €q )] dcp = iaW sin(a-e 0 ) 
a 

Regarding the result (20), it should "be noted that 
this formula naturally is not the only possible explicit 
presentation of the velocity function, since Cauchy's in- 
tegral formula can also be applied to every regular func- 
tion in the outside zone of the profile combined with 
i*(z). 

3. APPROXIMATE THEORY OP POTENTIAL PLOW 
OP THIN WING SECTIONS 



There is a theory for thin wing sections which, pro- 
ceeding from the image of a nonuniformly covered flat vor- 
tex layer, establishes an approximate relation between the 
velocity differences on the upper and lower surfaces of 
the wing section and the slope of its mean line (reference 
5). It constituted for a long time the only one affording 
a functional relation between velocity distribution along 
the wing section and its contour, and proved very useful 
despite its restricted range of . validity . Hence, the de- 
sire to expand this theory to include medium wing sections 
in approximate validity, though" no satisfactory solution 
has been found up to now..' 
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Theodorsen's theory ,.• described in the preceding sec- 
tion., affords the- strict?- relation ship between the wing- . 
section shape and the velocity distribution along the con- 
tour. It, moreover, affords information about thin" wing 
sections, and it might be conjectured that for vanishing 
thickn'ess' the equations of this theory with proper omis- 
sions, would become the equations of the old theory of 
thin airfoils, which would mean that Theodorsen's theory, 
represents the desired generalization of the old theory. 
However, that is not the case, as is readily proved: 
Theodorsen's theory contains a theory of thin wing sec-.- 
tions as limiting case, the results of which are a system 
of equations different from those' of the theory. of the . 
lifting vortex surface. • 

Jhe theory developed, hereinafter leads, for vanishing 
airfoil thickness, to the formulas of the linear .vortex, 
theory and represents, in its formulas, if not in thought, 
a generalization of the old theory. In distinction from 
Theodorsen's theory, it. is fundamentally an approximate 
theory but it involves in return no infinite method. Its 
result is rigorous only in the specific case of the flat 
plate exposed at angle a to the stream; but the more the 
wing section departs from the flat plate - that is, the 
greater its thickness and camber - the more the results 
assume an approximate . character . 

The method, like that of Theedorsen, begins with the 
cdnformal transformation of plane "z into a plane. z - 
with the given airfoil - by means of the- function 

, z = z + fe*~~± 

whereby the wing-section contour becomes a curve z (6 ) = 
\|/( G) + i6 

© , not much different from the unit circle. For 

this almost circular contour, that flow with the complex 

velocity w*(co). = W e~ 1(X at infinity must be determined, 
for which at point z = 1, which corresponds to the trail- 
ing edge, the velocity assumes the value zero. 

But, while in Theodorsen's method the flow about the 
almost circular curve is rigorously Computed by means of' 
the analytical function - obtainable, to be sure, only in 
infinitely many steps - which transforms the contour, of 
this curve smoothly into the curve of the unit circle,, the 
solution of this flow is now effected approximately oh 
the basis of the following reasoning: 
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At point P. : oh- : the contour of the almost circular 
curve the tangential" Velocity has a certain magnitude 
u t (6). Hence the' ; radial velocity component in P is, for 
small d approximately . u t ( 6) d(e) ? whereby the radial 
velocity in outward direction, and the- tangential velocity 
in direction of increasing 6 carry the positive sign 
(fig. 3 ).' 

About the- magnitude of the tangential velocity in P, 
the corresponding. flow ahout the unit circle affords, for 
the present, some information because, when 1 the almost 
circular curve changes into the unit circle in which dis- 
tance PQ, and angle £ uniformly, tend toward zero for 
all 8, the flow changes into that about the unit circle 
and converges, the tangential velocity . u t (G ) in P against 
the tangential velocity u tQ (6) at point' Q, of the unit 
circle . 

Then it is. assumed that the distance PQ, and the an- 
gle £ are so small for all .6 that no difference need 
be made between u' t (6) and u t (6) in the intended ap- 
proximate calculation. Hence the radial velocity in P is 

approximated at u+ ■& and, while the change on transi- 

o 

tion from P to Q. is ignored, this quantity substitutes 
for the normal velocity u n (e) at point Q, of the unit 
circle. Accordingly, since the flow about the unit circle 

(with velocity W e~* la at infinity and with the rear stag- 
nation uoint at TDoint z = l) is u t = - 2 W [sin(e-a) + 

" "0 

sin a] , we have : 

u n (G) =u t (e)*(9)=-2W [sin(e-a)' + sin a] *(6 ) 

* ° (21) 

In this fashion the problem approximately reduces 
to the second limiting problem of the potential theory for 
the unit circle-. Formulated for the complex velocity, the 
problem consists of giving a regular analytical function 
in the outside of the unit circle, which has the value 

¥ e~ la ' at infinity, a zero point at- z = 1, and whose 
normal component on the unit circle takes the value pre- 
scribed by equation (.21). Be sired , above all else,, is the 
explicit expression for the tangential velocity on the 
unit circle. 
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In any flow at point z = e*^ of the unit circle, the 
following equation's* ... . 

u n = u x cos G + tLy sin G, u t = - u^. sin 6 + u y cos 8 

exist "between the normal component u n (6), the tangential 
component u^ ( 6 ) , and the components u x (6) and u y (6) 
of the velocity in x and y direction. 

Combining u x and u y into the complex velocity 
function u*(z) = % - i Uy affords, for z = e 3 -^ 

[izu*(z)] = e i£ (i Ujj.+u y ) = (cos 6 +i sin 6)(i u x +u y ) = 

= - u x sin & +u y cos 6 + i (u x cos 6 +u y sin 6) u t ( G )+i u n ( 6 ) 

The result is a distribution of the tangential veloc- 
ity over the circle periphery feasible for a potential 
flow with the aid of 'Poisson's integral as real part to the 
imaginary part u n (G), according to equation (21). This 
function 

2" . 2rr 



JL_ 
2tt 
■o 



J u n (G') cot i--^- dG< = * j* [ain(6»-a) + 



G ' F 

+ sin a] (eO cot — -- — dG' 

is, however, not the desired result because, since the re- 
lated analytical function i z u*(z) at infinity is reg- 
ular, u*(co) =0 is contrary to the posed problem for the 
related flow. 



To obtain the desired flow, the most general poten- 
tial flow for which the normal velocity on the circle 
periphery disappears and the velocity at infinity has the 
correct value, must be superposed. This is the flow about 
the unit circle with the tangential velocity 

u t (6 ) = - 2w [sin(G-a) + sin 8] 

o *" 

Hence 
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u t (6) = - 2W |sin(G-a) + sin 0 - ~ J [sin(G'-a) + 



+ sin a] *(e») cot d6j 

Lastly, the available parameter 8 must he so deter- 
mined that the velocity disappears for 6 = 0 (Kutta- 
Joukowski theory). For 8 the equation reads: 

2tt 

sih p = sin a + ~ j [sin (6 '-a) + sin a] * (e i ) cot d6< 



which finally gives the sought-for .tangential velocity as 



u t 



( 6 ) = - 2W -f sin (6 -a) + sin a + 



2TT 



+ / [sin (6" -a) + sin a] 0 (8 ' ) 



(22) 



(cot —- - cot 



e » - e 



) * e *} 



or, transformed, 

sin (6' -a) + sin a = sin 8' (cos a + sin a tan --^ 



s m 



cot -x cot 



sin 



sin 9 » 



. e » .• • e '-e 

sm — g- sm — - — 



6 ' 

n ci COS ~ — , 

6 ' — — — — — — — — — 

sin 0 e S in 6 • . e «-e 

2 cos — sm — - — 



e > 

cos -p- 
sm_B 

6 6 ' - G 

cos - sm — 
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'and consequently 
■u t (6) = 2W .sin e 



cos a + sin a, tan — - 

2 



2TT 



2tt cos S ... 

2 o 




cosa+sina tan t^-^(e)' 



cos -- 



sin 



G '-6 



dG 



V(22a) 



This equation presents the velocity distribution on 
the almost circular contour in plane z; from it the dis- 
tribution over the airfoil contour is obtained by trans- 
formation of the z plane into the original plane with the 

aid oi the function z = ^ (z + ^ . The enlargement fac- 
tor of this transformation for z = e^ +i6 is 



dz 
dz 



* e""^ ^sin 2 $ + sinh a 



and the absolute velocity on the airfoil contour with W = 
2W (absolute velocity at infinity) is: 



|w(6)| = W 



e^ sin G 



J~ sin 3 G + sinh 2 \|/ 



2TT 

ri 



cos a + sin a tan — 

2 



COS 

(cosa+sinca tan— ) d (6 ' ) — dG 

* 2/ . fii-fi 



2tt cos 



6 J 
2 o 



sin — 



(23) 



In' order to establish a relationship between this re- 
sult and the linear vortex theory, the equation needs fur- 
ther conversion. To this end, the periodic function 



f(8) = ( 



cos a + sin a tan 



f)*(e) 



below the integral is divided in a component fi (6) symmet- 
rical to £ = tt , whose Fourier series consists of cos 
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terms, and an un symmetrical component f s (e) whose Four- 
ier series consists of sine terms only. If #1(6) and 
£ 2 (e) are the thus-declared components of function 6 ) . 
that is, if 

= ! [He )+*(2tt- e)], d a (6) =J [*(e )-*(2tt-6)] (24) 

we have : 

fi ( G ) = I [f (G)+f (2TT-6 )] = cos a « 1 (e)+Bin cx tan | * a (9 ) 

c 

fo (6) = i [f (G )-f (2-rr-e)] = cos a £ 2 (6)+sin a tan _ * x (6 ) 

2 ^ 

Owing to the symmetrical qualities of the functions 
f x (6) and f 2 (6 ) , the integral in equation (23) can "be 
substantially simplified. It is 3 



Auxiliary conversion formulas are: 

.' x+y . x-y 
x sin — i - sin -p- 
1 • 1 \ cos X 2_ 

2 I . x-y . x+y 2 . x-y . x+y 

\ sin — sin -~ J sin -g- sm -tj*- 



x 

cos 



2 co? ^ sin 2 2 cos 2 £ 

d cos 2 sm 2 y 2 

= r. O S — ■ =z er. sm — 

PI, \ 2 cos x - cos y 
- i (cos s - cos y) 



= - 2 sin ? 



1 + cos x 



cog £/ + _i_\ = cos J 



2 cos x - cos y 

sin 5£Z + sin 2£Z 



2 V -sin sin ~£ y sm sm -g- 

2- sin f cos ^ ,2 sin f cos § 

x. . 2 2_ 9 _ n „£. - - 

= cos - r— - . = -c;cos 2cosx~cosy 
- — (cos x - cos y ; 



y sin x 

= - 2 cos— 

2 cos x - cos y 
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e i e 1 

cos -- . ... P .. cos — 

2 j r ( / 4> f rs ,\ 2 



fl (e.) 4 d 6« = / f x (e») — ~t^"«' 

sin J S in 

d n - <- 



° 2tt - B 
r» cos k 



+ y fl ( 2w ^) — - 2 -~-r~e a(2n - e) = 

TT 2 

jr ., g I . e ' 

cos-- co s -p 

f. (m)/ ~ 1 d£ ' = 

sin - — a — — sin 



2 oj-" g 



TT 



6 / / \ 1 + cos 6' . c 
-2 sin ^ / ^ (CO ^sT~"cc7? 46 



o 

1 ikewise 



f 1 T 1 G i 

cos-- , ' cos — 

t ) * — ae > 

. G 1 - G 



cos — i 

^ (E ..__i_ 4£ fa(£ 



sin — — 0 sin 2 



o 



2tt - e 

co s — 



+ / f a (2 T T-"5') . ~: d(2TT-G) = 

J * sin 2tt - I -__£ 



TT 



2 



TT / 6 1 6 1 

TT /COS -S- COS ~ 

/ fa (e.)(-~— tt- e l a£! 

. \ sm - — _- sm g 



TT 
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( 



1 - X' 



3 2 

cosh \|/ 0 -x 



co see* since. 



•J 1 



~x\ 



+1 



1 +1 

TT 



«b" (| )d! 



+ 1 

ft 



-1 



co sa 



+1 



-1 



-i 



> 



u 1 v cosh \|/ u -x v 



cosa-sina 



l- x\ 
1+xJ 



+1 



+1 



p**&hL]-i cosa /!=£ /I±f jll 

! TT./ | - X J TT </ 1+X.J 1 V 1-| | -X 

-1 "I 



+ 1 

n 



+ — sina 

TT 



1+1 l-ac 



-1 



(26) 



This general result contains several specially note- 
worthy cases. 

1. Flat plate.- In the case of flat plate, it is 
v|/ 0 (x) eTy u (x). = . 0,. and also. 6 0 ( x ) = * u ( x ) = °> as a re ~ 
suit of- which equation (26) reduces : to 



[i(x)J = W 



cos oc. ± sxn a 



-_x 

+ X 



in accord with the previously derived equation (lO) for 
the absolute velocity on the plate surface. This is the 
only case where equations (26 ) give the rigorous result. 

2. Symme trica l air foil.- For this, it is ty Q (x ) = 
\|/ u (x) and 4 0 (x) = - * u (x); hence ^(x) EE 0, d a (x) = 
d 0 (x). In this case the velocity distribution is, accord- 
ing to equation (26): 



( 



|»(x)| =W / 
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.1 - x' 



cash 2 \J/ 0 ~x s 



cosa ± sma 




V !)<i£ 



(27) 



This equation is fundamentally an approximate equa- 
tion. First of all , . no good agreement can he expected for 
the trail ing-edge vicinity - that is, "for x values ap- 
proaching +1, if the edge angle is finite. In fact, for 
6 (l) ^ 0, the principal value of integral 



| - X 



increases "beyond all limits if x> 



1. 



The explanation 'for this is as follows: With finite 
angle at the trailing edge, the tangential velocity u^.(G) 
approaches zero differently for small 8 than it does at 
the circle. Hence it is not justifiable to solve the nor- 
mal velocity according to equation (21 ) on the oasis of 
the velocity distribution over the circle. And on this 
point, equation (27) must also fail. 

For practical purposes, it is recommended to figure 
either with a sharp trailing edge or, what amounts to the 
same, round off the. trailing edge and locate point z = 1 
midway between peak and center of peak curvature. In the 
second case, tHe' first factor of equation (27 ) , and with 
it the velocity, approaches zero for 6 = O". 

3i_Infini_tely_thin_airfo.il.- For infinitely thin air- 
foils, it is ^ 0 (x) = - \|/ u (x) and -b Q (x ) = d u (x); hence 

^i(x) = •6 0 (x), ^a(x) = 0, Since, in this case, the func- 
tion \|/ 0 (x ) disappears for x = + 1 and x = - 1, \J/ Q (x ) 
may be approximated at ^ c (x) 5 0 if the camber is slight, 

while function $ Q (x ) ~ J|(x) can, in this case, be ex- 



d6 

pressed with i> (x ) = - ^.Z - - y 

G ■• dx ■ 



because the airfoil 



follows completely in proximity of the piece of the x axis 
between -1 and +1. Accordingly, equation (26) gives: 



I 



2 8 
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|w(x)| =¥ 



cosa - — sina 

TT 



-1 



l=x 

1+x 



■ - 1 



sma- i cosa/ 

TT 



/Li _n 

V 1+ £ £~X 




1±L _d t 
1-1 S-x 



(28) 



The chief result of the linear vortex theory is a re- 
lationship "between the velocity differences on the top 
and "bottom of an infinitely thin airfoil and the slope 
y'(|) of its median line. Equation (28) gives for this 
relation the formula 

+ 1 



w 0 (x)-w u (x)=2f /IgLina-i cosa, J y'(£) /gf (29) 

-1 



which for a .= 0 reduces to 



+ 1 

n 



i.u^j.-a/s-y ,.( t ) 



-1 



X 



and represents the known result of the linear vortex the- 
ory (references*). The approximate theory developed here 
for medium thick airfoils contains the' old linear vortex^ 
theory as special case and, in that respect., represents its 
generalization, 

5. CONTRIBUTIONS TO GENERAL POTENTIAL THEORY 

OF.¥ING SECTIONS 



■ The potential flow about an airfoil is completely de- 
scribed "by the related complex velocity function, as ex- 
emplified in equations. (4), (9), and. (20), which give the 
velocity function flows about the unit., circle , the flat 
plate, and about any other airfoil.' The" velocity function 
is an analytical . function,, regular at. any point in the out- 
side zone of the' airfoil. For the vlcinijby of t he inf i- 

_ _. 



Fuchs-Hopf, Aerodynamik, Bd. 2, S. 87. 
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nitely remote point the series development (2a) is appli- 
cable. The coefficients of this series are associated 
with the force and moment applied by the flow on the body 
(equation 3). 

In the following, these arguments about the complex 
velocity are substantially complemented, wherein relation- 
ship is established between function and velocity distri- 
bution on the airfoil contour and with the direction of 
the airfoil contour. This reasoning leads first to devel- 
oping new integral equations for the basic problem of air- 
foil theory, which connect the velocity distribution along 
the airfoil contour with the direction of the airfoil con- 
tour at every point in exact manner; and secondly, it in- 
dicates a method for developing airfoils with required 
characteristics. 

Assume that the airfoil is plotted in the complex £ 
plane. If the flow strikes the airfoil at an angle a, 
then let w*(a; £) be the complex velocity. This function 
is regular in the outer region of the wing. For £ — ^ <* , 
it is : 

lim w*(a; O = W e~ ia 

The residuum of w*(a;£) is purely imaginary. For the 
values of w*(a;£) along the airfoil contour, the for- 
mula 

w*(a;£(s)) =.w*(a;s) =fw(a;s)| e ~^(a;s) (30) 

is valid, whereby s is the arc length on the contour 
measured from the trailing edge in positive direction of 
rotation )w(a;s)| indicates the absolute velocity at 
point s, and A(a;s) is the angle formed by the veloci- 
ty and the positive real axis (fig. 4). 

Further, let £ (z) represent the analytical function 
which transforms the contour of the airfoil in the £ 
plane into the contour of the unit circle in the z plane. 
The function £ (z) ' is regular for all finite „. z with 
|z| > 1. In the foregoing, l(z) had been standardized 

by the requirement that lim should be real, the im- 

z — >oo dz 

age point of the airfo il . trailing edge being .a wholly, ar- 
bitrary point of the unit circle. In distinction from it. 
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a subsequent rotation with the angle-- P ■ .. is. to bring the . 
image point of the trailing -edge to point z =1. Accord- 
ingly, it is at infinity ■ 

lim = as-iP 

Z >a> o-z 

For great z the series development 

/SB \ 
t(z)=a e- 1 ^ z(l+ — +-§ + ...) (P lt P a ••• complex) (31) 
V z z > 

Then the substitution £(z) in w*(a;£) gives an 
analytical function 

w*(a;z) = w*(a;£(z)) (32) 

which is regular outside of the unit, circle and assumes 
for z = e icp , according to equation (30 ), the values: 

[w*(a;z)] _ iC p = |w(oo;s(cp))| e- 1 * <<P > ) = 
z ■ — s 

= |w(a;cp)| e-^^T^) (30a) 

This function serves as "basis for the_subsequent investi- 
gations- The series development of w*(a;z) for great 
z is intended for later on. Because of equation (31) and 
of the development 

= C 0 + °£ + || + (o 0 - f e~^ C x = ££) 

valid for *w*(a;|-), the first terms of this series are: 

a z 



«iP / Pip 
w*(cg; z) = W e" la + ' C x 2— + ( C 2 -Cj. p x a e ) + . . . (32a) 

az \ f a z 



The behavior of the function w*(a;z) at the point 
z = + 1, resembling the trailing edge, is of special in- 
terest. If t 0 is. the point of. the trailing edge in 
plane £ and V the edge angle, we find: 

u - u ~ (z - i) TT - . or: ( Z . - i) ~ a - t 0 ) 
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for the vicinity of point z'= 1. 

Since in the flow about the circle with stagnation 
point at z = 1 , the complex velocity, acts like (z - l) 
the flow ah out the airfoil „ 

w*(a;0 ~ ( z - 1) U ~ t 0 ) 

For the function w*(a;z) resulting from suhstituting 
5 (z), it is: 

/ 2?t-y\ V Y 

i*(a;z) ~f( z ~ 1 ) 77 ) 2tt "Y _ ( z M i )" 

Hence the following result: If the airfoil has at the 
trailing edge the angle Y, then 

Y 

w*(a;z ) = const-. ( z - 1 ) TT (33) 

for the vicinity of z = 1. 

The relation of conformal function £(z) with 
w.*(a; £) • and w*(a;z) has some surprising • consequences . 
The function £(z) ties the complex -velocity of flow about 
the arbitrary airfoil to the velocity function of the flow 
about the unit circle. If W e~ icx ' is the velocity at 
infinity for the flow about the unit circle, and z = e icp i, 
the forward stagnation point, the complex velocity, since 
the rear stagnation point must lie in the image point of 
the trailing edge,, i.e., at z = + 1 .. is, according to 
equat ion (4 ) 

W e- ia ' (l - j).(l - --j~0 with q>i = tt + 2a' 
where, according to equation (6) 

Jt..,.-.-t-(l-l)(l.^)-^ , (34, 

From this it follows that the forward stagnation point 
of the flow about the unit circle is coincident with the 
zero point w*(a,;z) belonging to the forward stagnation 
point- of the airfoil flow. Herewith , cp x and, since cp x = 

tt + 2a 1 , a* itself is determined. for the infinitely 
remote point equation (34) gives; 
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— f JL ~ r (36) 
w* (a; z ; 



a e -3 ^ = W« e- ia ' W* 1 e ia 

hence 

a = |iy.3 = a'-a (35) 

The cenformal function £(z) itself follows from 
equation (34) "by integration. It 'is 

1 

which gives for z = e iCp the profile contour £(cp) re- 
lated to a predetermined function W*(a;z). : From the fact 
that the airfoil contour is a closed curve, it follows 
that the residuum of the integrand in equation (36) must 
"be equal to zero. The function w*(a;z) therefore, has 
the important characteristic 

U \\ r e i(TT+2a+2P)-] . 

Residuum Ul - ±) |l - & i I = 0 (37) 

I \ zj I z j w* (a; z ) J 

For the Coordination of the circle periphery and the 
airfoil contour there exists, according to equation (34), 
the following relationship: 



dcp 



sin 



(cp - a ' ) + sin a 



| w(a;cp)| 



(38) 



This equation is recommended for solving s(cp) when the 
distribution of the absolute values of w*(a;z), that is, 
the function |w (a;cp)| on the unit circle is known. a' 
must "be so chosen that I w (a;rr+2a 1 ) I = 0. Then the func- 
tion 1w(a;cp)l can "be replotted- with s(cp) and the ve- 
locity distribution |w(a;s)| on the profile contour ob- 
tained. 

Conversely, for given velocity distribution |w(a;s)| 
equation (38) give's, inasmuch as point s x of the forward 
stagnation point oh the profile contour corresponds with 
point cp = tt + 2a' 

s cp 

^ |w(a;s)| ds = 2¥ «| j |sin(cp-a') + sin a'J dcp or 

s -i Tr+2a t ■ 

s 

j |w(a;s)| ds = 2W'| co s (cp-a, ' ) + cO s a '+ (TT+2a '-cp) sin a '| (39) 
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wherein the constants <x ,: and V are left for determi- 
nation. This is achieved- >y" the fact' that s = 0 relates 
to cp = 0 and s = S to cp = 2tt, where S is the air- 
foil contour. The, solution for a' and W is: 



s. 



|w(cc;s)| ds 1 + (?± + a ») tan a 



n(a') 



|w(a;s)| ds 1 - tan a* 



I w ( a ; s ) | d s 



4 cos a ' + 



(I 



+ a 



sin a 



t 



► (39a) 



The function M-(a') 
creases "between a' 



occurring in the first equation in- 
a' = + 2 monotoni cally 



- ?f and 



from 0 to op. Hence there is one, and only one a' for 
each value of the quotient "between these limits. Then the 



second equation gives the value 



(39) gives with these values 
s(cp), with, the aid of which 
"be obtained from |w(a;s)| . 



a ' 

the 



W for 

and W ' 
f unct ion 



a 1 . . Equat ion 
the function 
|w(a;cp)| can 



Hence the remarkable fact that for the conformal trans- 
formation "between the contours of an arbitrary airfoil and 
the unit circle, the coordination of the "boundaries can "be 
determined without knowing the shape of the profile if the 
velocity distribution along the profile contour for an ar- 
bitrary potential flow about the profile is known. 

Lastly, the distribution |w(a;s)| for any other an- 
gle can be deduced according to equation (38), equally 
without knowing the profile, from the velocity distribu- 
tion |w(a;s)| for a certain angle of attack a. Since 
the same function s(cp) relates to both flows about the 
same profile, it is with cp = cp(s) 
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|w(a;s)| = |w(a;s)| 



sin(cp~a '-a+a) + sin(a'+a-a) 



sin (cp-a 1 ) + sin a f 



(40) 



In particular, it gives for a = a — a 1 = - 3 the veloc- 
ity distribution for the circulation free flow 



w ( s ) | = | w (a; s )| 



sin cp 



sin(cp-ot') + sin a' 



(40a) 



Thus the velocity distribution along the profile con- 
tour related to a certain potential flow about the profile 
determines the velocity distribution for other. angles of 
attack as well, and can be obtained by elementary calcula- 
tion. And the fact that the related profile itself can be 
defined in simple manner without having recourse to an in- 
finite method is one of the ensuing results. 

Following these more general arguments about the ve- 
locity function and the velocity distributions on profile 
contour and circle periphery, the relationship with func- 
tion d(s) and '•6(cp) decisive for the shape of the con- 
tour, is discussed. Reverting back to equation (30a), we 
replace cp by s, so that it now reads: 

[w*(a;z)] iC p = |w(a;cp)- |- e -l*(a;<p) (30a) 

What are the connections between the functions |w(a;cp)/ 
and $ (cc;cp ) ? 

With the complex function •'w*(a;z) new complex func- 
tions f(z) = f ( z ,w* (a; z ) ) can be formed which are regu- 
lar in the outside zone of the unit circle with inclusion 
of the infinitely remote- point, and on the circle contour. 
With' ' 

f ( z , w* (a; z ) ) = g(u,cp) + i h(u,cp) 

the two . interrelated equations exist between real and im- 
aginary parts of f(z) on the unit circle 
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g(i,cp) = e(»).. 



h(l,qp) = h(») + 



at the same time^ g (l , cp) ..and h(l,cp.) are functions of 
cp^ w(a;cp) and d(a;cp), according to equation (30a). The 
result therefore is a diversity of integral equations, all 
of which establish the relationship "between the functions 
lw(a;cp)| and $(a;cp) in exact manner. A few illustra- 
tions are given : 

1. The function w*(a;z) has itself the characteris- 
tics required of f(z). The next example therefore is 
f(z) = w*(a;z). In this case 

g(l,cp) = | w | cos g(o») = W cos a 

and 

h(l,cp) = - |w| sin d, h (on.) = - W sin a 

Hence the identities 



sin "5 (op') cot ^ ^ ^ dcp' 



cos •?(<¥>') cot -2 dcp" 

Every such equation can fundamentally "be used to de- 
termine the other function for given |w(cp)J or 1(cp), 



2tt 



27T , h('l ,(p'y cot 5L-~^ dcp' 



2TT " 

y g (i f cp«) cot seljp-ff dcp« 

" o • 



|w(cp)l cos #(cp) = W cos a + 

2TT 

i- f is(»'>i 

0 

|w(<p)| sin d(cp) = W sin a - 

2-TT 



2. ./ 
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and it should "be possible in various ways to give infinite 
iteration methods which lead to results, and' might even 
De advantageotis for the application. This finding like- 
wise explains the surprising fact incident to the compar- 
ison of" the two theories (3 and. 4) that different assump- 
tions may lead' to substantially unlike profile theories 
which in limiting cases, as of the infinitely thin air- 
foil, give outwardly entirely different equations. 

2. A very interesting example is the following case: 
f(z) = | (z - |) (w*(a;z) - w*(co) 
Here Qf (z)] lcp = i sin cp(w e - J * - W. e~ ia ) =. 

= Iwl sin !i sin cp - W. sin a sin. Cp.+ 
+ i(|w! cos ¥ sin cp. ~. W cos a sin cp) 

ir 

and, on account of equation 32a) with Ci = g— 



= iCel- - _E_ (- sin 3 + i cos p) 
4-rra 4-Tra 

Hence the two identities read: 

|w| cos i sin cp — W cos Ou sin cp — 

2tt 

Si 

r COS 6 1 / / i — I . T • i 

= 4na + 2n J ( 1 wl Sin * Sin * * ~ 

o 

- W sin a sin cp ' ) cot dcp 1 = 



- W sin a cos cp + 

n 

+ ~ I |w| sin sin cp 1 cot 3?~ r_JP dcp' 
o 



£(») = 



1 lm g 
z — &co 



r_cos__§ 

4-rra 
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The value of the constant L~ c -0.s_.£ follows from the 

4iTa 

Kutta- Joukowski condition that at the trailing edge - that 
is, for cp = o, the' velocity must he finite. It is 




o 



hence the integral equation: 



| w C cp) J cos #(cp) = W Coos ot + sin a tan + 



2tt 



(41) 




o 



This exact equation is intimately related to the the- 
ory of the lifting vortex surface. For the replacement of 
lw(cp)| below the integral by the velocity distribution of 
the- fla-t plate, according to equation (10 ) in first approx- 
imation, leads to_an equation which (substituting 6 for 
cp, and d for ■&) agrees with equation (23) for the case 
\|/ = 0 . Thus the result is again among others, the result of 
the theory of the lifting vortex surface - this time, how- 
ever, not as limiting case of a somewhat more generalized 
approximate theory as in (4), but as first approximation 
of an exact equation for the most general flow about any 
arbitrary profile. 

3. The equations obtained, according to the described 
principle are, as a rule, integral equations possessing, 
as in_the two examples, combinations of functions |w(cp)| 
and •S(cp) below and outside of the integral. However, 
theso functions can be separated by so choosing the func- 
tion f(z) that its real part depends upon |w(cp)| , and 
its imaginary part on #(cp) only. 

If the airfoil has the trail ing-edge angle Y , _ and 
the forward stagnation point of the flow lies at e 3 -^ 1 , 
then . according to equation (33), 



f (z,w*(ot;z)) 




is a regular function with the characteristic required for 
all J z | S 1 . 
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Now in general 



In 



/ i A \ 

(i - V) 



_ 2 = ei<P 



= In 



(2 sin * 



for A "g cp < 2tt + X, while for cp < X and for cp > 2tt+ X 
the two functions for real and imaginary part must "be con- 
tinued with the period 2tt. Hence (for X = 0 and X = 



In 

L 

and 



d ~ l)] = In (2 sin %) + i Hrffi for 0 < cp < 2tt 

^ z >J z = e ic P V ^ • 



r t e iCPl X" 
In ( 1 - ~ 



icp 
z=e ^ 



.' |«Px - <*> . 
= In 2 sin ) + 1 



- 0 5r cp < cp x 



•) + 1 c— — + 2) 



for •( 

v cp x < cp < 2tt 



In this case it gives for f(z) = g + i h 



g(l,q>) = In 



I 2 sin I J (2 sin — 4j j 



g(oo)=ln IT 



(l)C p) = .^(cp)^ 3£2-*£?± I (for cp<^, h(co) = - 



a 



hence Poisson's integral ' leads to 



|w(cp)l =W (2 sin I) 



(2 sin - 
\ 2 



2 IT 

o / ■ 
|-9x-CP|_^ 2tt 0 L 



cot M dcp' 



(43) 
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= ■ a - - —J- 2 - < £±-f-% ± 'I - 




Jzl£P_!iLl_ cot" d^' ■ (44) 



( 2 S in f /(a sin i^Pl) 



2 



The advantage of using the function In w*(a;z) as 
"basis had already "been recognized "by F. Weinig, who like- 
wise cited equations (43) and (44). 

The discussion of. the general airfoil theories is con- 
cluded with Weinig^ In w method, involving the question 
of profiles of prescribed characteristics, of profiles re- 
lating to a given velocity distribution, and the potential 
flow about a given profile. 

5. THE WEIJTIG- In w METHOD 3 



•In preparation of the treatment of profiles with pre- 
determined characteristics, the appearance of the most im- 
portant profile characteristics with the function w*(a;z) 
is first ascertained. For this purpose equation (42) is 
used to formulate w* (a; z ) with F(z) = f(z) - f (°° ) 

Y 

w*(a;z) = We" 1 * (l - l)" (l - e F < z > 



with 



F(z) = Z = £ ^nli^a, since F(co) = 0 

n = i z n n*=i z 11 



(45) 



The parameter Y in the formula is, according to equa- 
tion (33) the trailing-edge angle of the airfoil. The fact 
that the contpur is closed affords some prediction concern- 
ing the coefficient c x . If it relates to an airfoil with 



2 _ — — — 

The connections between the velocity distributions at dif- 
ferent angles of attack and the In w method has been re- 
ported by F. Weinig in the following publications: Z.f.a.M.M. 
9 (1929), p. 507; Werft Reed. Hafen 14 (1933), p. 131; Luftf. 
Forschg. 12 (1935), p. 221. 
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fixed center of pressure, a second condition for c 3 must 
"be added. The condition for' c x follows from equation 
(37). It is 

Residuum (l - (l - ^7 t) = 

V z/\ z / w* (a; z )J 



= e.™ r e s i duum * 



l- 7 - -Z ^ 



/ . i\ * n=l z n e la / . . Y \ 

0-i) e |= -w-C- 1+ n~ c 0 = 



hence the final condition 



a, = - 1 + T)j. = 0 (46) 

The condition of fixed center of pressure follows 
from the fact that an airfoil has a- fixed center of pres- 
sure when the moment for zero lift is zero. In this case 
the aerodynamic center of v. Mises' lift paraoola is zero, 
and the air force therefore passes through the aerodynamic 
center of the airfoil at every angle of attack (reference 
5). This means, according to equation (3), that "by series 
development for w*(cx;£) in the case of C x = 0 for 
C 2 = A 2 + i B 2 , the relation 

A 0 B 3 + B 0 A 2 = 0 

must "oe fulfilled. 

According to equation (32a) the disappearance of Gi 
is accompanied "by that of the first coefficient- C x of 
the series development for w*(a;z), and the second coef- 
ficient G 2 is, in this case, connected to C 2 "by 

' C 2 = a 2 e- 2i P C 2 . 

But , according to equation (45 ) 

00 / ■ Y e inc *M 1 
In w*(a;z),= In W - i a + (^c n - — —J ~ E 

hence , 



w*(cc;z) = W e~ ia |l + (c x -'J - e'i«Pi) \ 

♦[(•» -5;- + -I- el '0° ^ + -J} 



+ 

2 
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That is, 

V'-'* e - ia (c-'-'J - e 1Cp i) = - W e" 3 - 00 (1 + e iqP M 

therefore "Cj. = 0! if cp x = tt; that is, if a = - 8. This 

is the condition for lift-free flow about the wing sec- 
tion . 

For e*^ 1 = - 1 and a = - B, we have: 

C 3 - T e ** (c 2 - - §) - W e±P (a 2 + i h 2 - ^ - J) 

C a = a - W e-*P (a a + i b a - "g. - l) 



hence 



The condition for fixed center of pressure is A 0 B 2 + B 0 A 3 = 
0, which means that the imaginary part of C 0 C a must dis- 
appear. But, since a = - 3, the constant' is C 0 = W e 1 ^ 

C o 0 2 = a 2 (a 2 + i b 2 - — - |) 

so that the condition for airfoils with fixed center of 
pressure reads 

b 2 = 0 (47) 

The closing condition and the condition of fixed 
center of pressure are also predictions about the func- 
tions w(a;cp) and •S(a;cp). The particular formulas can 
be readily obtained on the identity: 

[In w*Ca; z).] z=e icp = In |w(a;<p)| - i £(a;cp)" = 

'= ln'W - i a + E ( c n _ _1 - S^lY-incp- 

n=i \ nTr n ✓ 

by division into real and imaginary parts and integration. 
It is for function In w(a;cp) ' '" 



2tt 

rt 



In |w(a;cp)| dcp '= 2tt In W. (48a) 



42 



■Ttf.A.C.A. Technical . Memorandum. Ho . -886 



2tt 



2tt 



In |w(a;cp) sin cp dcp = ~ tt sin cp x 



r (closing (46a) 
condi.t i on ) 



-/ In | w(a;cp)| cos-cp dcp = - tt(1+cos ep x ) 



w(a;cp) sin 2cp dcp = - g sin 2cp x 



and for function ^(a;cp) 
2-rr 

^ $(a; cp) dcp' = 2tt a 
0 2tt . 



£(a;cp) sin cp dcp = - tt(1+cos cp x ) 



2tt 



J £(a;cp) cos cp dcp = + tt sin cp x 



2tt 



TT 



•6(a;cp) cos 2cp dcp = + - sin 2cp x 



(fixed cen- 
ter of 
pre ssure ) 



( clo sing 
condition ) 



(fixed cen- 
ter of 
tires stire ) 



(47a) 



(481) ) 



(461) ) 



(471) ) 



The general theory presented here is superior to the 
old method, part icularly" in the manner of- developing math- 
ematical profiles. The calculation needs no special con- 
formal functions with a number of correctly chosen parame- 
ters as oasis, hut rather proceeds from the general- func- 
tion w*(a;z) which is connected with the airfoil charac- 
teristics- in the described manner. 

Either the coefficients c n in equation (45 ) can ho 
prescribed, noting that Cj must always = - 1 + -, and 



that t) 2 = 
f unct ion 
tack a, 



0 for fixed center of pressure; or else, the 
|w(a;cp)| can he given for a certain angle of at- 
from which the velocity distributions along the 
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profile contour (and hence the flow characteristics, of the 
airf oil in question ) ean he deduced, according to equation; 
(38) and (40). By the selection of |i(a;cp)| , it should 
he observed that equations (46a) are complied with. and, if 
a fixed centor_of pressure is involved, equation (47a). 
For cp = 0 , |w(cx,;cp)| must,, like ep' y/!rr , become zero if 
the airfoil is to have the t rail ing-edge angle Y. 

With this chosen function lw(a; 



In 



lw (a ; cp)| 



Y 



(a sin |7 ( 8 sin M) 

must be formulated and developed in Fourier series. If 
this development reads 

In LSisiffiiL = in + 

(2 sin ff (Z sin M) 

CO 

+ Z (a. cos a cp + k sin n cp) = 
n=i n n 

2tt 

In w(a;cp) dcp - '^1 - — J cos cp + 



2tt . y 
o 



+ £ (a n cos 'ii cp + h n sin n cp) 



n=2 



the related function ^(ajcp) is, according to equation 
(44): 



+ £ (a n sin n cp - Td_ cos n cp) 
n =2 

and the profile contour itself, according to equation (36), 
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- 1) .0 



i^A dz 
" z J w*(ct; z) 



with z = e ic P. Hereby, according to equation (4), 



cp 



z = e ic P. 



,icp 



y . . . dz = y • . . i dcp, 
1 © 



[w*(a; z )] z = Q icp = |w(a;cp)| e 



-I*Ca;q>) , 



[,...-»'.(x-l) (i-^) 



z-e ic P 



= 2i W» e~ iCp (sin (cp-a') + sin a') 
and consequently, ^ 

U<P) = £ - 2a W / ll^S£=SlL±-^i2L_SLi e i^(a;^) dc P (36a) 
0 "" -/ jw (a; cp)] 

o 

The separation of real and imaginary parts gives for the 
nrofile contour the nararaeters 

£(cp) = £ 0 -2a W / linCgg^pj^o* cos ^ (a;cp) dcp 
</ |w(a;cp;| 
o cp 



Tl(cp) = Tl 0 -2a W J 



8i£(fe^aO+l^_SL! S in ^(a;cp) dcp 
| w (co;cp)| 



(36b) 



The constant a is a scale factor for the profile. These 
integrals can he graphically evaluated. 

As regards the question of profile for a' given velocity 
distribution jw(a;s)| along the profile contour, raised 
previously, the following can now he stated. If it has 
previously "been proved that |w(a;s)| belongs to a poten- 
tial flow about a closed profile, then |w(a;cp)| can be 
ascertained from equation (39), and the cited method gives 
the profile contour. Hence the problem of determining the 
profile for a certain velocity distribution is, as origi- 
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nally discussed by Weinig, explicitly solvable and conse- 
quently essentially simpler than the reversed problem. 
Nevertheless, the function jw(a; s )| is no suitable start- 
ing point for profile calculations, - for the choice of 
lw(a;cp)| as the choice of |w(.a;s)|* is restricted by the 
final condition (46a). These equations are, as follows by 
substitution of cp ( s ) , a very complex requirement for 
function |w(a;s)| which is. not summarily met, and would 
involve a return to |w(a;cp)| . which would bring us back 
to the previously described method. 

Lastly, the problem of computing the velocity distri-_ 
bution for a prescribed airfoil can be treated by the In w 
method. The deciding condition is equation (43). But, 
since with given profile (and forward stagnation point) the 
function 6(a;s) is exactly known, while d(cx;cp) is not 
known, the problem leads progressively to the executed ap- 
proximation of the function $(<x;cp) which, however, means 
that the principal problem of the airfoil theory leads 
here also to an infinite method. 

The first formula- for £(a;cp) can be obtained by the 
following process: 

Select the junction line between trailing edge and 
airfoil nose as profile axis. Draw a straight line from 
point P(s) on the contour to this line. If the base 
point P« has the coordinate x, and t is the profile 
chord, put 

/ \ x(s) 
cp(s) = arc cos ^J£~ 

and coordinate the upper surface of the profile to the in- 
terval 0 < cp < tt , and the lower surface to interval 
tt < cp < 2tt (fig,. 5). 

How, £(a,;s) is replotted on the basis of this approx- 
imately valid coordination of s and cp; let, thereby, the 

function £(ct;cp)_ originate. Since s(cp) is not accurate- 
ly known, this #(a;cp) generally does not comply with the 
final condition (46b), and therefore is not as yet practi- 
cal for the calculation, according to equation (43). For 
this reason, a correction of the f orm • A-b (a;cp ) = A a sin cp 
- Ab_cos cp with the free constants Aa and Ab , is added 
to £(a;cp) in such a way that 
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&i (a;cp) ' = '£(ajcp) + A £(oc.;cp) 

complies with equation (46b). ¥ith this l 1 (a;cp) a first 
approximation. I Wj. (a;cp )| of the velocity distribution is 
then attempted by means of equation (43). 

This correction in respect to (46"b) resultsin the 
suppression of the cos term in the Fourier series for 

a(a;cp) + - E-— - + 5Pl_T_SS + 2 and the addition of the co- 
tt 2 2 2 

efficient - 1 + - to the sine term. Then the series 

TT 

for the integrand of (43) reads: 

*x(a S cp> + J + + g = - (l - J) sin cp + 

+ E (a n sin n cp - b n co s n cp ; 
n=2 

and the first approximation |w 1 (a;cp)| is, according to 
oquat ion (43 ) : 

■ 1 

TT 

| w\ (a;cp)| = W'(2 sin |) - x 

~( 1 - - X ;coscp + E (a„ cos ncp +b n sin n cp ) 
X (2 sin -^-> ^ ^ n=S " 

With this function f w x ( <x ; cp ) i the connection between 
s and cp can then he newly computed, according to equa- 
tion (38), and a new approximation ■6 2 (a;cp) obtained from 

•6(a;s). The method can be continued at will. Its con- 
vergence is not easily followed mathematically. But there 
is good reason to surmise that results are quickly achieved, 
and also confirmed by Weinig. 

Translation by J. Vanier, 
National Advisory Committee 
for Aeronautics. 
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Figs. 1,2,3 
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Figure 2. 




Figure 3. 
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Figs. 4,5 





Figure 5. 



